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Abstract. 

An attractive way to generate neutrino masses as required to account for current 
neutrino oscillation data involves the spontaneous breaking of lepton number. The 
resulting majoron may pick up a mass due to gravity. If its mass lies in the kilovolt 
scale, the majoron can play the role of late-decaying Dark Matter (LDDM) , decaying 
mainly to neutrinos. In general the majoron has also a sub-dominant decay to two 
photons leading to a mono-energetic emission line which can be used as a test of the 
LDDM scenario. We compare expected photon emission rates with observations in 
order to obtain model independent restrictions on the relevant parameters. We also 
illustrate the resulting sensitivities within an explicit seesaw realisation, where the 
majoron couples to photons due to the presence of a Higgs triplet. 
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1. Introduction 

While solar and atmospheric neutrino experiments [U [2, [3] are confirmed by recent data 
from reactors [I] and accelerators j5] indicating unambiguously that neutrinos oscillate 
and have mass [6], current limits on the absolute neutrino mass scale, 

mv < 1 eV (1) 

that follow from beta [7] and double beta decay studies [8], together with cosmological 
observations of the cosmic microwave background (CMB) and large scale structure [9] 
preclude neutrinos from playing a direct role as dark matter. 

However, the mechanism of neutrino mass generation may provide the clue to the 
origin and nature of dark matter. The point is that it is not unlikely that neutrinos get 
their mass through spontaneous breaking of ungauged lepton number [101 [11]. In this 
case one expects that, due to non-perturbative quantum gravity effects that explicitly 
break global symmetries [12], the associated pseudoscalar Nambu-Goldstone boson - 
the majoron J - will pick up a mass, which we assume to be at the kilovolt scale |13j . 
The gauge singlet majorons resulting from the associated spontaneous L-violation will 
decay, with a very small decay rate T, mainly to neutrinos. However, the smallness of 
neutrino masses (Eq. (CD)) implies that its couplings to neutrinos gj vv are rather tiny and 
hence its mean life is extremely long, typically longer than the age of the Universe. As 
a result such majorons can provide a substantial fraction, possibly all, of the observed 
cosmological dark matter. 

Here we show how the late-decaying majoron dark matter (LDDM) scenario, and 
in particular the majoron couplings gj vv and gj 77 to neutrino and photons respectively, 
can be constrained by cosmological and astrophysical observations. 

The paper is organized as follows. In Sec. [2] we describe the basic cosmological 
constraints on the LDDM scenario, in Sec. [3] we describe the "indirect detection" of 
the LDDM scenario and determine the restrictions on the relevant parameters that 
follow from the x-ray background and the emission from dark matter dominated regions. 
In Sec. H] we compare the sensitivities of CMB and x-ray observations to the LDDM 
scenario, stressing the importance of the parameter R = Tj 77 /r j vv . Finally in Sec.[5]we 
discuss an explicit seesaw model realisation of the LDDM scenario, where the majoron 
couples to photons due to the presence of a Higgs triplet. 

2. Cosmological constraints 

The LDDM hypothesis can be probed through the study of the CMB anisotropy 
spectrum. In fact, current observations, mainly of the Wilkinson Microwave Anisotropy 
Probe (WMAP) lead to important restrictions. Indeed, the LDDM scenario has been 
explored in detail within a modified ACDM cosmological model; in particular, it has 
been shown that the CMB anisotropies can be used to constrain the lifetime tj ~ Tj^ v 
and the present density Qj = pj/p c of the majoron [15]. 
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The reason is that the late decay of majorons to neutrinos would produce too much 
power at large scales, through the late integrated Sachs- Wolfe effect, thus spoiling the 
CMB anisotropy spectrum. WMAP third year data [Hj can be used to constrain: 

r Jvv < 1.3 x 10~ 19 sec _1 , (2) 

at 95% CL p2]. 

This result is independent of the exact value of the decaying dark matter particle 
mass, and is quite general, in the sense that a similar bound applies to all invisible 
decays of cold or warm dark matter particles [T6l \T7\ [T8] . 

The CMB spectrum can also be used to constrain the majoron energy density. This 
can be translated to a limit on the majoron mass in a model-dependent way. Given the 
majoron mass raj and lifetime tj, the present majoron density parameter Qj can be 
written as: 

^ 2 = ■ 3 T£tv^' TJ < < 3 > 

where h is the dimensionless Hubble constant, to is the present age of the Universe, and 
the parameter (3 encodes our ignorance about the number density of majorons. The 
normalization in Eq. [3] is chosen such that f3 — 1 if (i) the majoron was in thermal 
equilibrium in the early Universe; and (ii), it decoupled sufficiently early, when all the 
quantum degrees of freedom in the standard model of fundamental interactions were 
excited. 

This simple picture can be changed if: (i) The majoron could not thermalize before 
it decoupled from the other species, or (ii) The entropy generated by the annihilation 
of some particle beyond the standard model diluted the majoron abundance after its 
decoupling. 

In any case it is reasonable to assume that the majoron decoupled at T > 170 GeV 
since its couplings to all the other particles in the standard model (SM) are tiny. Using 
the WMAP third year data, the following constraint on Qjh 2 can be obtained (95% 
CL.) assuming the dark matter to consist only of majorons [15j : 

0.09 < Qjh 2 < 0.13 (4) 

Since Eq. fl2]) implies tj ^> t , the above constraint together with Eq. ([3]) gives: 

0.12 keV < (3m j < 0.17keV (5) 

Our ignorance of the details of the majoron production mechanism, namely the 
value of (3, can always be used in order to accommodate additional restrictions to the 
majoron mass mj coming from observations of the large scale structures. 

The limits quoted above apply to the invisible decay J — > vv. There exist also 
the very interesting possibility to use the CMB polarization to directly constrain the 
radiative decay J — > 77. This is because photons produced by dark matter decay 
can inject energy into the baryonic gas and thus affect its ionization history. This 
will ultimately lead to modifications of the CMB temperature-polarization (TE) cross- 
correlation and polarization auto-correlation (EE) power spectra. In Ref. [19], WMAP 
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third year data are used to obtain the following constraint for the radiative decay width 
r ra d of long lived dark matter particles like the majoron: 

Cr rad < 2.4 x lO-^sec" 1 , (6) 

where £ is an "efficiency" factor describing the fraction of the decay energy actually 
deposited in the baryon gas. This depends, among other things, on the energy of the 
emitted photon. As a rule of thumb, consider that for redshifts 10 < z < 1000, when 
most of the hydrogen is neutral, photons with energy in the range 13.6 eV (the hydrogen 
ionization threshold) to approximately 1 keV will transfer most of their energy to the 
baryon gas through photo-ionizations, so ( ~ 1. On the other hand, the Universe is 
transparent with respect to the propagation of photons with E ^> 1 keV, and in this 
energy range one expects to have ( ~ and then no significant upper limit on r ra d can 
be obtained from CMB polarization. 



3. X-ray analysis 

In a variety of neutrino mass generation models with spontaneous violation of lepton 
number, majorons have an effective interaction term with photons 

(jj'yyJ 6 ^ Fy^Fpu . (7) 

Majorons in the keV range are therefore expected also to decay radiatively into two 
photons of energy ~ mj/2, since the decay can be considered to a very good 
approximation as happening in the dark matter rest frame. This leads to a mono- 
energetic emission line as a characteristic signal of our decaying dark matter model. 

Such an emission line could be possibly be detected both in the diffuse x-ray 
background and in the emission from dark matter dominated regions. We now consider 
the constraints coming from both kinds of observations. 

In the following, when necessary, we will consider a LDDM scenario within a ACDM 
cosmology with Q\ = 0.75, QoAih 2 = 0.11, Qbh 2 = 0.022, and h = 0.72, corresponding 
to the best fit values of the CMB analysis in Ref. |15j . 



3.1. Diffuse x-ray background 

Photons produced in late majoron decays will show up in the diffuse x-ray background, 
if the Universe is transparent with respect to their propagation. This is indeed the 
case after the Universe has been completely reionized (z < 10): photo-ionization is no 
more effective in absorbing the photon energy, simply because there are no more neutral 
hydrogen atoms to be ionized. 

The flux F(E) of decay radiation at the present time (z = 0) is given by [20] : 

c / E \ 3 N^T j^nj(z) 



F(E) 



4tt \E 7 J H(z) 



(8) 

l+z=E 1 /E 



where A 7 = 2 is the number of photons produced in each decay, nj(z) is the number 
density of majorons at redshift z, H(z) is the Hubble parameter, and E 1 = mj/2 is 
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the energy of the photons produced in the decay. This can differ from the energy E 
at which we are observing due to the cosmological redshift of photons. In other words, 
when looking today at an energy E < mj/2 we can still expect some signal from photons 
emitted in the past with energy E y = nij/2 that have been red-shifted to lower energies. 

However we know from the CMB that the majoron is very long lived, so that we 
expect the decay spectrum to be dominated by very recent decays. We model the 
spectrum as mono-energetic with E = E 1 and a flux given by: 

F(E ) - ° N ~* Tj ^ n °' J (a) 
W ~4tt Ho ' [) 

where the subscript denotes quantities evaluated at the present time. 

This should be compared with the observed diffuse x-ray flux from ASCA [2T] and 

HEAO-1 [22], operating in the 0.4-7 keV and 3-500 keV ranges respectively. The flux 

can be modeled as [23] (units are sec -1 cm -2 sr _1 ): 



Fobs(E) 



E * ~ - 4 



keV 



1 keV < E < 25 keV, 



E x - 1 ' 6 



380 — , 25 keV < E < 350 keV, (10) 
keV , 



E ^ -°- 7 



350 keV < E < 500 keV, 



Below 1 keV, the strong galactic emission must be carefully removed in order to 
find the extragalactic signal [24] , and consequently we do not extrapolate the above 
approximation to lower energy for the purpose of the present analysis. 

Then, requiring E(E^) < E b s yields an upper limit for the majoron decay width 
to two photons T j 77 . In particular, in the range, IkeV < E 1 < 25keV we have: 

sec" 1 ~ v°- 72 y v O- 11 / ^ keV ' 

This limit, together with the constraints at higher energies, is shown in Fig. [TJ 

This simple analysis, and the resulting constraint, can be improved in two ways. 
First of all, one can look for small distortions in the smooth diffuse flux produced by a 
DM emission line that is possibly lying well below the background signal. In addition, 
one can take into account the contribution to the signal coming from the Milky Way. 
This was applied to the HEAO-1 data in Refs. [251 [26]; i n this way, the above constraints 
can be improved by as much as three orders of magnitude (see below). Finally, we note 
that bounds in the soft x-ray region can be obtained from the observations of a high- 
resolution spectrometer |27j . 



3.2. X-rays from dark matter dominated regions 

Observations of the x-ray emission from dark matter dominated regions can be used to 
restrict the decay rate into photons and the mass of any dark matter candidate with 
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a radiative two-body decay. This follows from the consideration that the detected flux 
from a dark matter dominated object gives a very conservative upper limit on the flux 
generated by dark matter decays in that object. 

Since the dark matter in cosmological structures is practically at rest (v/c w 10 -4 ), 
the line broadening due to motion of the dark matter is negligible compared to the 
instrumental resolution of current day x-ray detectors. Hence, a good instrumental 
spectral resolution increases the sensitivity to a mono-energetic emission line. 

For majorons the 0.1-0.3 keV x-ray interval is very interesting. Unfortunately this 
range is not accessible with any of the standard CCD instruments on board the present 
x-ray observatories Chandra and XMM-Newton. However, Chandra carries the High 
Resolution Camera (HRC) which combined with the Low Energy Transmission Grating 
(LETG) makes it possible to obtain spectra in the 0.07-10.0 keV range. The resolution 
of grating spectra is very high {Efwhm ~ 1 eV) [28] but all spatial information about 
the photon is lost, except that it is known to origin from within the field of view (apart 
from minor effects of scattering along the line of sight). 

To reach the maximum resolution requires bright point like sources located at 
the aim point of the observation. Unfortunately dark matter structures have spatially 
extended distributions and will produce a very faint signal, if any. Extended sources 
can be thought of as made up of many point sources, but then most of the sources 
are off-axis. The effect of a source being off-axis is, that in the detector plane, there 
is an ambiguity between angle and photon energy, which gives a "smearing" towards 
lower energies and hence a line broadening in the obtained spectra. The line smearing 
is energy dependent and worst for high energies [29|. 

No optimal sources for a search for dark matter decay line emission have been 
observed with HRC/LETG. Still, from grating observations of an active galaxy, we have 
improved the upper limit on the decay rate from the dark matter halo in which the 
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active galactic nuclei is embedded by orders of magnitude. For photon energies above 
0.3 keV better constraints are obtained from conventional x-ray CCD observations of 
merging clusters of galaxies such as the Bullet Cluster [30] or Abell 520 [29] . 

We have studied a Chandra HRC/LETG observation of the Seyfert 1 galaxy NGC 
3227 (observation id 1591), shown in Fig. El NGC3227 has a redshift of z = 0.004, 
which corresponds to a luminosity distance of 16.7 Mpc. The data has been processed 
and analysed with CIAO version 3.4 using CALDB version 3.4.0 [31 J . The obtained 
spectrum is shown in Fig. [2j 
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Figure 2. The observed Chandra HRC/LETG spectrum of NGC3227 (folded with 
the instrumental response). 



In general, the instrumental response cannot be unfolded from the spectrum in a 
model independent way. Instead a model is folded with the instrumental response and 
fitted to the data using x 2 statistics (here we have used the spectral fitting package 
Sherpa distributed with CIAO). 

The model is used to determine an upper limit on the received flux. Since no 
physical quantities are derived from the empirical model, it is chosen to fit the data 
(and as such do not necessarily represent a physical model of the emission). The data 
were split into two intervals: 0.072-0.276 keV and 0.276-4.14 keV and fitted separately 
to models composed of a power law and four Gaussians for the lower interval and two 
power laws and two Gaussians for the higher interval. In order to ensure that no emission 
lines are sticking above the model, the fitted model was re-normalised so there were no 
bins in the spectrum at more than 2a above the model. 

As mentioned above, any emission line is smeared out because of the instrumental 
resolution. Towards lower energies the smearing has the shape of a Gaussian with 
the width given by the instrumental resolution. Towards higher energies, where the 
extension of the source plays a role, the smearing depends on the overall distribution 
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of the dark matter. We have assumed an Navarro, Frenk & White (NFW) [32] profile 
for the dark matter halo of NGC3227 with a scale radius of 15 kpc and a virial radius 
of 200 kpc. These are conservative representative values for galaxies, in the sense that 
most galaxies have a smaller scale radius and a smaller radius, leading to less smearing 
(and thereby tighter constraints). The resolution is only sensitive to the full width at 
half maximum of the density profile, so that choosing a different parameterization (e.g. 
Moore [33]) would not significantly alter the results. 

Since we are only interested in the upper limit on the measured flux, it has been 
determined in slices of width £ 7 + FW HM instrumental > E > - FWHM smearing 
instead of the exact shape of the smeared lines (the difference between the two methods 
is negligible [29]). 

The mass of NGC3227 has been taken to be 10 11 solar masses which is conservatively 
low based on the luminosity of the galaxy [M]. The observational field of view is ~ 25 kpc 
at the distance of NGC3227, reducing the observed mass to about a tenth of the total 
mass. This is probably an underestimate of the observed mass, but a larger observed 
mass will only improve the constraints. 

Assuming only one kind of dark matter, the observed flux, F b s , at a given photon 
energy yields an upper limit on the decay rate from two-body radiatively decaying dark 
matter: 

r Jry < L . 12 

J77 - M fov 1 ; 

The determined flux is dominated by the baryonic emission of the galaxy, which varies 
with energy, introducing an apparent energy dependence on the constraint. 

The resulting constraint on the decay rate is shown in Fig. [3] together with earlier 
published constraints. 



4. X-ray versus CMB 

In the previous sections we have shown how the CMB can be used to constrain the 
invisible decay J — > uu, while x-ray observations can constrain the radiative decay 
J — >• 77 [j. From a theoretical point of view, a very important quantity is the branching 
ratio of the decay into photons BR(J — > 77) that, as long as the decay to neutrinos is 
by far the dominant channel, is given by the ratio of the decay widths: 

B R(J - 77 ) ~ R = £in. (13) 

J- Jvv 

Note that the decay J — » vv arises at the tree-level, while the radiative J —> 77 
mode proceeds only through a calculable loop diagram. Beyond this, theory can not 
predict the expected value of R, which is strongly mo del- dependent. Here we use R as a 

I We do not consider the CMB polarization limit in the following because (i) it depends on the efficiency 
of the energy transfer to the baryonic gas and (ii) it turns out to be less constraining than x-rays in 
the regions of interest. 
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Figure 3. Upper limit on the decay rate from NGC3227 (red), the Milky Way halo 
observed with a prototype cryogenic spectrometer (salmon) [27] , XMM observations of 
the Milky Way (sand) [35] . Chandra observations of the Bullet Cluster [3D] and M31 
[3JEU2] (orange), HEAO-1 observations of the diffuse x-ray background (aquamarine) 
[HI HSJ, INTEGRAL SPI line search in the Milky Way halo (blue) [Ml EH- Filled 
regions are excluded. 



phenomenological parameter varying over the wide range 10 -25 — 10~ 3 (see, for instance, 
Fig. [6] below). It should be clear that the observations described in Sees. [2] and [3] restrict 
Tj uu and Tj 77 , leaving the branching ratio unconstrained. 

However, we are also interested in knowing for which models the x-ray observations 
can probe the decaying majoron dark matter hypothesis with higher sensitivity than the 
CMB. In particular, we expect that models with large branching ratios will be better 
constrained by the x-ray observations, since they will predict a larger production of 
photons. 

The x-ray limits presented in Sec. [3] have a mass dependence, which we need to 
take into account in our assessment of the relative constraining power of the two types 
of observations. We know, however, from the CMB (see Eq. [5]) that: 



0.12 keV 0.17keV 



(14) 



(3 ~ ~ (3 

at 95% C.L. Fixing the value of f3 is then equivalent to fixing the majoron mass, apart 
from a small uncertainty (which we take into account, see below). We express our result 
in terms of j3 instead than mj. 



In order to compare the CMB and the x-ray constraints, we fix the value of (3 and 
determine the corresponding observational ratio of Tj 11 /Tj l/1 y. According to the CMB 
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Figure 4. Sensitivity of CMB and x-ray observations to the LDDM majoron scenario 
as a function of (5 and R = Tj 71 /Tj u ^. The black lines are the loci of points where 
R = R* , i.e., where the CMB and x-ray constraints (from a given object) are equivalent. 
In the hatched regions above the lines, X-ray constraints are stronger; below, CMB 
constraints are stronger. The color codes are the same as in Fig. [3] See Sec. 0] for 
discussion. 



constraints, we take the mass of the majoron to be equal to mj = 0.145 keV//3, with an 
associated la error of oj = 0.01keV//3. Then, we find the maximum Tj 77 allowed by 
the x-ray emission for this value of the mass (as explained in Sec. 13. 2p . The uncertainty 
in the exact value of the mass is taken into account by convolving the upper limits 
shown in Fig. [3] with a Gaussian of mean mj and variance equal to o"j. Let us call this 
value Tj^.We also denote with Tf^ = 1.3 x lO^sec" 1 the CMB upper limit on the 
decay width to neutrinos. Then for the following value of the branching ratio: 

pmax 

R* = -^L (15) 

1 Jvv 

the two sets of observations yield exactly the same constraining power. In other words, 
for this particular value of the branching ratio, it would be the same to constrain the 
decay rate to photons using the x-rays and then obtain the decay rate to neutrinos 
using Tj uu = Tj-y-y/R, or to do the contrary, i.e. to use the CMB to constrain the 
invisible neutrino decay channel and from that obtain a bound on the photon decay. 
Larger branching ratios (R > R*) will be better constrained by observations of the x- 
ray emission, while smaller branching ratios (R < R*) will be better constrained by the 
CMB. 

We repeated above procedure for f3 ranging from 10~ 5 to 1, comparing the x-ray 
constraints of Fig. [3], one at time, with the CMB constraint. We did not include the 
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XMM observations of the Milky Way because they are discontinous and this makes 
the mass-averaging procedure problematic. The results are illustrated in Fig. HI We 
can roughly say that for small majoron masses {(3 ~ 1), we should resort to x-ray 
observations to probe the region R > 10~ 6 , while we should use the CMB for R < 10~ 6 . 
For large neutrino masses, x-ray observations are better when R > 10~ 8 , while CMB is 
more informative for R < 10~ 10 . 



5. Particle physics 

We now turn to the particle physics of our decaying dark matter scenario. Although 
many attractive options are open [40J possibly the most popular scheme for generating 
neutrino masses is the seesaw. The simplest type I seesaw model has no induced majoron 
radiative decays. For this reason we consider the full seesaw model, which contains a 
Higgs boson triplet coupling to the lepton doublets [H]. 

In addition to the SM fields one has three electroweak gauge singlet right-handed 
neutrinos, i/£., a complex SU(2)l scalar triplet A, with hypercharge 1 and lepton number 
—2 , and a scalar singlet a, with lepton number 2. We will denote the scalar 577(2) 
doublet as 0. The Yukawa Lagrangian is given by 

Cy = r„Ql<K + Y d Ql<f>*d c L + Y e L T L <f)*e c L 

+ Y v L T L cj)vl + Y L L T L AL L + ^fu c L u c L a + H.c. (16) 

In order to extract the relevant couplings of the majoron that are responsible for the 
decays in Eq. (j2J), we review here the main steps of the procedure developed in [TT] . 
using the basic two-component Weyl description of neutrinos as in [JT]. 

Using the invariance of the scalar potential under the hypercharge U (l)y and lepton 
number U(1)l symmetries and assuming that these are broken spontaneously by the 
vacuum configuration, one finds, from Noether's theorem, the full structure of the mass 
matrix of the imaginary neutral component of the scalars given in terms of their vacuum 
expectation values (vevs) as [11] 



M 12 = C 



( 4 4 _ 2 ^L 2 a \ 

vi V2.V2 V2 

2 2 

2 v i v i vi_ 

2 ui v± 



;i7) 



\ fa v 3 / 

with C = d 2 V jd^o 1 and v\, v 2 , v 3 are the vevs of the singlet, the doublet, and the triplet, 
respectively. One sees that M 12 has a non zero eigenvalue, m 2 A = Tr M 12 and two null 
eigenvalues. These correspond to the Goldstone bosons eaten by the Z gauge boson 
and, as expected, to the physical Nambu-Goldstone boson associated with the breaking 
of U(1)l, the majoron J. The parameters of the scalar potential of the model can be 
chosen so that the pattern of vevs obtained by minimization respects the so-called (type 
II) seesaw form, namely [11] 

^3 < v 2 <^ Vi 



X-ray photons from late- decaying majoron dark matter 



12 



In particular, since the smallness of the triplet vev arises through the vev seesaw relation 
one can show that it is not spoiled by one-loop radiative corrections. 

The resulting profile of the majoron, J, following from Eq. ( TlTl) takes a very simple 
form in this seesaw approximation, namely [11] 

J^_M O/ + ^3 A O/ + a O/_ (lg) 
ViV 2 Vi 

In the presence of the gravitationally induced terms that give mass to the majoron, 
the mass matrix of Eq. ( TlTl) is slightly modified, but these effects are sub-leading and 
negligible. 

We are now ready to determine the coupling of the majoron with the light neutrinos. 
From Eq. ffl~6]) one obtains the full neutrino Majorana mass matrix as 

M u — - ( l LV3 ) . (19) 

2 ^ » 2 Y RVl J 1 ; 

so that the effective light neutrino Majorana mass matrix is given by [TT] : 

(20) 



The coupling gj vu of the majoron to the neutrinos can also be obtained using Noether's 
theorem according to the procedure described in Ref. [TT]. In this way one finds that 
the majoron couples to the mass eigenstate neutrinos proportionally to their masses, 

*. = -^- PD 

where V\ describes the scale at which the global lepton number symmetry breaks, 
typically 10 6 - 10 9 GeV (see below). 
The decay width T Juu is given by 

mj £ r (r<) 2 



Jvv 



(22) 



32tt 4v 2 ' 

Let's now turn to gj 77 . From the Yukawa Lagrangian of Eq. ( fi~6i) and from Eq. (JT 
we have that the majoron interacts with the charged fermions through 

2v 2 - 2v 2 

- — -Y f (-2T 3f )fj 5 f J = — ^-m f (-2T 3f )fj 5 f J , (23) 

V 2 V\ V 2 V\ ' 

where T 3 f is the weak isospin and we have assumed that the charged fermion mass 
matrices are diagonal. The interaction term of Eq. (123]) gives rise to the interaction 
term with photons given in Eq. ([7]), with an effective coupling given by 

ol „ , r-w-\ \ o I 1 m 2 T \ 2v 2 , , 

gjll = —X f N f -2T 3f Q} 1 + --i . 24 

2tx \ ^ m f } v z Vl 

where one notices the cancellation of the "anomalous- like" contribution S f Nf (—2T 3 f) Q 2 . 
As a result we have 
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with 

£ = 1 "fa 

7 E f N f (-2T 3f )Qj±^ vi ' 

where Q f and Ay are the electric charge of / and its colour factor, respectively. 

Fig. [5] shows how the currently allowed range of neutrino masses selects an allowed 
strip in the plane V\ — mj consistent with neutrino oscillation data [6] and with the 
cosmological bounds on neutrino mass [12], assuming that the CMB bound (T5]) on the 
J — ► vv decay rate is saturated. 




mj (keV) 

Figure 5. The strip indicates the region in the v\ — mj plane allowed by current 
neutrino oscillation [S] and cosmological data [42j , assuming the maximal J — > vv decay 
rate. The vertical lines delimit the mass values required by the CMB observations, for 
different values of j3. 

The lower lines correspond to the cases of normal and inverse hierarchical neutrino 
masses, while the top line holds when the three neutrinos are (quasi)-degenerate. The 
vertical bands in the figure indicate the mass region of Eq. [5] singled out by the CMB 
observations, for two different values of (5. 

We also note from Eq. (1251) that, for a fixed value of the majoron mass and the 
lepton number symmetry breaking scale v%, the two-photon decay rate only depends on 
the vev of the triplet and on the sum of the squared masses of the neutrinos, namely 
on the two possible scenarios in the neutrino sector, hierarchical or degenerate. For a 
given scenario the decay is then fixed only by v 3 , as it can been seen in Fig. [6], where 
the top panel corresponds to the hierarchical case while the bottom one holds for the 
quasi-degenerate spectrum. The diagonal lines in Fig. [6] give the dependence of T j 77 on 
mj, for different values of v^. As it can also be seen from Eq. ([251) . the largest values 
of correspond to the largest radiative decay rates. 
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mj (keV) mj (keV) 



Figure 6. Majoron decay rate to photons as a function of the majoron mass mj, 
for different values of the triplet vev, V3. We assume the invisible decay bound 
to be saturated. The top and bottom panels refer to hierarchical and degenerate 
neutrino mass spectra, respectively. The shaded regions are excluded by observations 
as described in Sec. [3l The vertical lines are the same as in Fig. [H 



One sees that in both scenarios small mj and v% values lead to decay rates well 
below the observational bounds. However, for large values of V3, say, v 3 = 5 GeV, 
roughly corresponding to the maximum compatible with precision measurements of 
electroweak parameters [13] , the radiative rates fall within the sensitivities of the Milky 
Way observations displayed in Fig. [31 and would be thereby observationally excluded. 
For lower masses the observational sensitivities would need to be improved by about 20 
orders of magnitude requiring completely new techniques from what is available today. 
The small radiative majoron decay rates would be avoided in models where the anomaly 
does not cancel due to the presence of extra fermions. We mention also in this case the 
possibility of further enhancement due to cumulative effects as those that might arise, 
for example, in higher dimensions. 

6. Summary 

We have investigated the production of x-ray photons in the late-decaying dark matter 
scenario, and quantified the sensitivity of current observations to such a mono-energetic 
emission line. In particular, we have studied the constraints from the diffuse x-ray 
observations, as well as by considering the fluxes generated by dark matter dominated 
objects. These observations provide a probe of radiative dark matter decays and can be 
used as an "indirect detection" of the LDDM majoron scenario. 

We have illustrated this explicitly for the case where neutrinos get mass a la seesaw, 
where the majoron couples to photons through its Higgs triplet admixture. Alternative 
particle physics realizations of the LDDM scenario can be envisaged, an issue which will 
be taken up elsewhere. Let us also mention that Majoron dark matter decays can be 
possibly probed in the future through 21-cm observations (see Ref. [H] for an application 
to other DM candidates). 
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